Study of the Gamma-ray Spectrum from the Galactic Center in view of Multi-TeV 

Dark Matter Candidates 



Alexander V. Belikov 
Institut d'Astrophysique de Paris, UMR 7095 CNRS, 
Universite Pierre et Marie Curie, 98 bis Boulevard Arago, Paris 75014, France 



Or 



Gabrijela Zaharijas 

Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151, Trieste, Italy and 
Institut d'Astrophysique de Paris, UMR 7095 CNRS, 
Universite Pierre et Marie Curie, 98 bis Boulevard Arago, Paris 75014, France 

Joseph Silk 

Institut d'Astrophysique de Paris, UMR 7095 CNRS, 
Universite Pierre et Marie Curie, 98 bis Boulevard Arago, Paris 75014, France 
Beecroft Institute of Particle Astrophysics and Cosmology, 
1 Keble Road, University of Oxford, Oxford 0X1 3RH UK and 
Department of Physics and Astronomy, 3701 San Martin Drive, 
The Johns Hopkins University, Baltimore MD 21218, USA 
(Dated: July 11, 2012) 

Motivated by the complex gamma-ray spectrum of the Galactic Center source now measured over 
five decades in energy, we revisit the issue of the role of dark matter annihilations in this interesting 
region. We reassess whether the emission measured by the HESS collaboration could be a signature 
of dark matter annihilation, and we use the Fermi LAT spectrum to model the emission from SgrA*, 
using power-law spectral fits. We find that good fits are achieved by a power law with an index 
~ 2.5 — 2.6, in combination with a spectrum similar to the one observed from pulsar population and 
with a spectrum from a > 10 TeV DM annihilating to a mixture of bb and harder t + t~ channels 
and with boost factors of the order of a hundred. Alternatively, we also consider the combination 
of a log-parabola fit with the DM contribution. Finally, as both the spectrum of gamma rays from 
the Galactic Center and the spectrum of cosmic ray electrons exhibit a cutoff at TeV energies, we 
study the dark matter fits to both data-sets. Constraining the spectral shape of the purported dark 
matter signal provides a robust way of comparing data. We find a marginal overlap only between 
the 99.999% C.L. regions in parameter space. 

PACS numbers: 95.35.+d 95.85.Pw, 98.70.Rz 



I. INTRODUCTION 

In the last few years, a number of new generation 
cosmic-ray detectors have provided high quality measure- 
ments of gamma rays (Fermi LAT [H, VERITAS 0, 
HESS H, MAGIC ||), charged cosmic rays (PAMELA 
0, CREAM @, Fermi LAT, HESS) and neutrinos (ICE 
CUBE 0, Antares @) at GcV and/or TeV energies. The 
newly acquired data are currently revolutionizing our un- 
derstanding of the origins and propagation of Galactic 
cosmic rays and they are equally important for indirect 
searches for dark matter (DM) . The theory of cold dark 
matter is supported by numerous pieces of evidence (e.g. 
an impressive agreement between predictions of N-body 
simulations and the observed large-scale structure of the 
Universe, Big-Bang Nucleosynthesis, gravitational lens- 
ing). Weakly Interacting Massive Particles (WIMPs) are 
dark matter candidates with the required cosmological 
properties and are naturally present in many beyond- 
the-Standard-Model theories, as a new particle of mass 
close to the Electro Weak (EW) breaking scale. An ad- 
ditional appeal of this paradigm is the so-called WIMP 
miracle: a TeV particle weakly coupled to the Standard 



Model (SM) particles that decouples thermally in the 
Early Universe, naturally has an abundance compara- 
ble to the DM abundance measured today. If it were to 
decay or annihilate to Standard Model particles, such as 
photons and electrons/positrons, the bulk of the emitted 
energy is expected to fall just below the DM mass, and 
to show features in this energy range. 

The trajectory of gamma rays, unlike that of charged 
cosmic rays, is unaffected by the Galaxy, allowing the 
analysis of individual sources. In order to maximize the 
luminosity of gamma rays produced in DM annihilations, 
a standard search strategy is to look at nearby sources 
with high DM density. An example of such 'good tar- 
get' is the center of our Galaxy. It is one of the closest 
targets and it is expected to contain a high dark matter 
density, due to a possible cusp in the density profile pre- 
dicted by N-body simulations [9HJl|. However, the dy- 
namical center of our Galaxy hosts a supermassive black 
hole (compact radio source SgrA*) [12], whose spectrum 
is measured over a range of wavelengths, spanning radio 
and sub-m m flp near-infrared 0, Ell, X rays fl6l [Iff], 
and 7 rays and which presents a complex back- 

ground for DM searches [25[ . 
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The Very High Energy (VHE) gamma-ray source at 
the Galactic Center (GC), HESS J1745-290, was detected 
by the HESS collaboration in 2003 and observed since 
[ll [22|, in the range 200 GeV through 55 TeV (an early 
detection in VHE gamma rays was done by the WHIP- 
PLE [H and CANGAROO-II telescopes). The GC 
is also observed by MAGIC [H] and VERITAS 0] with 
spectra compatible to the one measured by HESS. The 
large zenith angles for these experiments however, result 
in significantly smaller exposures. The gamma ray emis- 
sion from the GC was subsequently interpreted both in 
terms of standard astrophysical emission originating in 
the vicinity of the supermassive black hole |26| - |3G | as in 
terms of a signal from self-annihilating dark matter 

The Fermi satellite, launched in 2008, measured the 
emission from the GC in the lower energy range 100 
MeV-100 GeV. In the first year Fermi LAT catalog 1FGL 
the source closest to the GC, was classified as 1FGL 
J1745.6-2900c. While the most recent data from this 
complicated region have not yet been analyzed by the 
collaboration (however see Fermi LAT proceedings and 
talks (36l-l38j), a dedicated analysis has been done by sev- 
eral authors [39h42| . In these works it was found that 
the Fermi and HESS sources are spatially coincident (in 
agreement with the preliminary findings of the Fermi 
LAT team) , though the issue of the extension of the GeV 
source is still unsettled. 

At energies above 10 GeV, the Fermi LAT has a good 
angular resolution comparable to that of HESS at TeV 
energies. This motivates the joint analysis of the two 
data-sets which, taken together, show the possible pres- 
ence of one or more spectral breaks. In [43J, it was 



shown that the spectra is well modeled by two plateaux 
with a steeply falling region in between, and astrophysi- 
cal mechanisms which rely on modeling collisions of en- 
ergetic protons, accelerated by the supermassive black 
hole, with the surrounding gas have been proposed in or- 
der to match this spectral feature (for an earlier model 
which accommodates a break between the GeV and TeV 
photons originating in the vicinity of the SgrA* see [2(| ) . 

Motivated by the complex gamma ray spectrum of the 
Galactic Centre source we revisit the issue of a multi- 
TeV dark matter fit to this region. We show that a sin- 
gle power law over the GeV and TeV range, with two 
bump-like features, provides a comparably good fit. We 
assume the feature at high-energy as a signature of DM 
annihilation. We then notice that the lower energy bump 
is curiously located around ~ 3 GeV, an average energy 
at which a harder emission from the pulsar population 
has a high energy cut-off 0, EH • In addition, we also 
show that a log-parabola in the low energy range (typ- 
ical for spectra of some AGNs (4f|) provides the same 
conclusions. 

In other words, we go back to the question of whether 
the emission measured by the HESS collaboration could 
explained by DM annihilation, but we use the Fermi LAT 
spectrum to model the standard astrophysical emission 



from the SgrA* 1 . 

It has been known from early works (3"H - [33T | that high 
energy HESS data, if interpreted as DM signatures, could 
be due to > 10 TeV DM, which is somewhat heavier than 
the standard expectation for a WIMP, and that the flat- 
ness of the HESS spectrum over two decades in energy is 
hard to fit with the usual WIMP spectrum. The second 
issue is somewhat alleviated in this work by the addition 
of a power law constrained by the Fermi LAT data. 

In this work, we focus on spectral features of the cen- 
tral source and leave morphological studies aside (for re- 
cent work addressing this issue see [H, E§|), which is ex- 
pected to be a g ood approximation for strongly cusped 
DM profiles [5C)t [5l| . While there are indications that 
the GC source at GeV energies can be spatially extended 
[4(il - |4"2| , this issue is still unsettled due to the complexity 
of this active region and we do not address it here 2 . 

On the side of spectral signatures of heavy ( > 100 
GeV) dark matter candidates, it has been recently noted 
that addition of the emission of W or Z bosons from final 
and internal states, neglected in standard c omp utations, 
could substantially alter the spectrum [53l - l58j . both at 
low energies and in the energy range close to the DM 
mass, which is critical for the spectral fits in the GC 
region. As the magnitude and shape of the spectral cor- 
rections in the hard energy region are model dependent 
we do not consider them here. We however comment on 
the implications of such corrections on our fits in section 

EDS 

Recently there have been strong indications of a line 
signature in the Fermi LAT data at around 130 GeV, 
with a likely origin within < 4° within the GC l59l - [6l| . 
The morphology of this emission was studied in [HJ find- 
ing a best fit position ~ 1.5° off the Galactic Center, 
with a ~ 5tr significance of the signal compared to the 
null-hypothesis. As the nature of this signature is still 
not clear at the moment, we do not address it here. We 
simply notice that if the signal proves to be of a non- 
instrumental origin a DM interpretation would be well 
motivated (but see [H, [64| for a difficulty to reach needed 
branching ratios to line channels in the usual WIMP 
models and [65| for alternative astrophysical explanation 
of such signal) . 

In addition to the gamma-ray data analysis, we make 
a connection with another astrophysical measurement 
which recently has shown a surprising feature: the lo- 



1 For work along similar lines, which instead uses the HESS data 
to model the standard astrophysical signal in the GeV range see 

mm. 

2 Note that at TeV energies a more extended region, ~ 2 deg 
(200 pc) around the GC, known as the Galactic Centre Ridge 
1531 also has a puzzling spectrum, harder than the surrounding 
medium. This emission correlates spatially with a complex of 
giant molecular clouds in the central 200 pc of the Milky Way 
and one possible explanation, based on energetics of the emission, 
is that it could have come from a single supernova explosion 
injecting hard cosmic rays, around 10 4 years ago. 
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cal population of electrons/positrons as measured by 
PAMELA, Fermi LAT and HESS experiments (see the 
recent review in (6(| for both dark matter annihilation 
and astrophysical interpretations). While it has been 
shown already that dark matter is viable but not nec- 
essary to explain these data sets, we return here to the 
DM interpretation in an attempt to do a more careful 
calculation for heavy candidates. We consider the full 
rclativistic calculation of energy losses as implemented in 
the GALPROP code USUI, critical for such high ener- 
gies, and, as opposed to the gamma-ray case, we account 
for electroweak bremsstrahlung contributions to the elec- 
tron and positron spectra, which can affect positron fits 
performed at low energies. Naively, as the high energy 
cut-off in the HESS gamma-ray data is higher than that 
measured by the HESS electron data, these data imply a 
different mass of the potential DM candidate. We nev- 
ertheless explore this possibility by deriving the best-fit 
DM regions, marginalizing over DM and astrophysical 
modeling. 

The article is organized as follows: in the section |TT] 
we briefly summarize observations on the GC source fo- 
cusing on gamma-ray spectra; in section fill Al we review 
DM density profiles in the central parts of our halo and 
annihilation spectra of heavy DM candidates in MI Bl In 
section UVl we describe the fits we perform to the Galactic 
center gamma-ray data, while in section [V] we focus on 
the electron and positron measurement. We compare our 
findings with the current constraints in the literature on 
multi-TeV DM models in I VII and summarize in section 

Em 



II. ASTROPHYSICS AT THE GALACTIC 
CENTER 

Below we briefly summarize current observations on 
the GC, and some of the proposed emission models in 
the gamma-ray range. 

Emission from Sgr A* in the radio to sub-mm presents 
a very hard luminosity spectrum L u ~ v a with an index 
a = 6.8 with a cut-off at about v ~ 103 GHz [nj]. The 
angular size of the source depends on the frequency of 
observation: at 1 GHz, it is of the order of 1.5 arcsec 
while it reduces to 0.2 milliarcsec at 86 GHz. Variabil- 
ity (though without a clear pattern) is observed at these 
frequencies. 

The quiescent flux in the near-infrared has been de- 
tected with the VLT [6{|, in the energy range 0.3-1 eV, 
as a point source with a position coincident with the su- 
pcrmassivc black hole within an accuracy of 10-20 marc- 
ses. Chandra X-ray observatory [j^], covering the energy 
range 0.1-10 keV with an angular resolution of 0.5 arcsec, 
observed the GC as an extended source with a size of 1.5 
arcsec. The near-infrared and X-ray emissions from Sgr 
A* arc characterized by a large variability (on different 
timcscales in the two cases). In addition to the central 
source, continuum X-ray and radio observations indicates 



the central 10 pc region of the Galaxy (comparable to the 
resolution of gamma-ray detectors) as a complex environ- 
ment characterized by several distinct structures [H, [71| . 

In VHE gamma rays, the latest analysis of the 93 
hours of live time observation of the central source by 
the HESS collaboration [22j revealed that the spectra of 
HESS J1745-290 cannot be fitted with a single power 
law, but that instead a power law spectrum with an 
exponential cut-off characterized by a photon index of 
2.10 ± 0.04 (stat) ± 0.10 (syst) and a cut-off energy 
at 15.7 ± 3.4 (stat) ± 2.5 (syst) TeV, provides a good 
fit. A dedicated analysis in [72| was able to localize the 
position of the TeV point source within 7.3 arcsec ±8.7 
arcsec (stat) ± 8.5 arcsec (syst) from Sgr A*, consistent 
as well with a pulsar wind nebula 8.7 arcsec from the 
GC [Z3[ 3 . Unlike in the radio and X-Ray observations 
of the GC, no variability has been observed in the high 
energy regime [22| . potentially indicating that the 7-ray 
emission mechanism differs from the low energy regime. 

The authors of [43[ analyzed 25 months of the Fermi 
LAT data using the P6V3 instrument response function, 
modeling the diffuse emission and point sources in the 
10° x 10° region and in the 100 MeV-300 GeV energy 
range. The spectral curve for the central source 1FGL 
J1745. 6-2900 was derived to be a power law with a slope 
r = 2.212 ± 0.005 and a flux normalization F = (1.39 ± 
0.02) x 10~ 8 cm" 2 s" 1 MeV" 1 at 100 McV. The authors 
noted that a better fit is obtained when the spectrum is 
described by two distinct power laws. In this case, the 
fitted slope is equal to T = 2.196±0.001 in 300 MeV-5 
GeV energy range, and T= 2.681±0.003 in 5-100 GeV 
energy range . It was also shown in this work that Fermi 
LAT source 1FGL J1745.6-2900 lies within the error box 
of HESS source J1745-290. 

Various models predict VHE emission from the cen- 
tral source, produced either close to the black hole itself 
[26l [30| or within an ~ 10 pc zone around Sgr A* due 
to the interaction of run-away protons with the ambient 
medium ^$2^. 

A theoretical model proposed in [43| which can explain 
the shape of the combined Fermi LAT and the HESS 
spectrum considers the hadronic interactions of relativis- 
tic protons which, having diffused away from a central 
source, presumably the central black hole, fill the inner 
few parsecs of our Galaxy. The size of the region com- 
bined with the energy dependence of the diffusion then 
result in a specific spectral shape. Authors show that the 
spectra can be well fit after several fitting parameters 



3 An angular resolution of approximately 0.1 deg corresponds the 
distance of 15 pc, when Solar position at 8.5 kpc is assumed. 

4 It has also been noted in |43H that the error of the lowest data 
point might be underestimated. Fermi LAT has a very broad 
PSF at lower energies, starting from 4° at 100 MeV and rapidly 
improving to 0.1° at > 10 GeV. Thus, taking into account the 
possible source confusion in the region, one should treat the first 
point in the spectrum (100 - 300 MeV) with caution. 
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are adjusted (free parameters in this model character- 
ize geometry, diffusion and injection rate history). This 
model has the additional appeal of reconciling variabil- 
ity of SgrA* at lower energies [74[ with no variability in 
gamma rays (as the gamma ray emission originates from 
a larger region surrounding the black hole and is emitted 
during the diffusion of the rclativistic protons through 
the interstellar medium). However, due to our poor un- 
derstanding of the spectrum of protons accelerated by 
the central black hole, or of past activity of SgrA*, it is 
not simple to independently test the best-fit parameters 
in this model. 

In this work we remain agnostic about the astrophysi- 
cal mechanism responsible for the emission of SgrA* and 
we sample different fitting functions to the Fermi LAT 
data in order to explore the uncertainty it brings into the 
DM interpretation of the HESS measurement . 

III. DM SIGNAL: DISTRIBUTION AND 
SPECTRA 

A. Signal 

The flux of gamma rays from annihilating dark matter 
in the direction ip can be written in the following way: 




where M x is the mass of the WIMP, (av) is the standard 
thermal cross-section (assumed to 3 x 10 -26 unless oth- 
erwise stated), p(r) is the dark matter density as a func- 
tion of distance r = \/l 2 + D 2 — 21D cos ip away from the 
center. The generalized NFW density profile is given by 
p(r) = (r/a )T( 1 + T . /a) 3-T ! where a is a length scale, p is 
the density scale and 7 sets the inner slope. The total 
flux for a source observed in a solid angle Af2 is 

<f> x (E) = ( dfl^ x (E^) 
J An 

d£7 7 v J 8irM 2 J A J lo / V ; 
= <t> x (£)J(AQ). (2) 

The particle physics and energy dependence is contained 
in the function $> X (E), while the astrophysics is contained 
in the factor J. We assume the spectrum dN/dE 1 to be 
given by a sum over annihilation channels 

dE 1 ^ dEi' [ ' 

We assume that AO = 10~ 5 sr (corresponding to the 
definition of a point source given the 0.1° resolution of 
HESS). As we perform only spectral fits, the details of the 
profile shape are not critical. In what follows we assume 
the standard NFW dark matter density profile @ with 



the following normalization parameters: the Milky Way 
mass Mmw = 1-6 x 10 12 M Q [73, the dark matter density 
Pq(Rq = 8.5 kpc) = 0.3GeVcm -3 and inner slope 7 = 
1.0. 

In the following we will express part of our results in 
terms of a 'boost factor' Bp. We define it as an enhance- 
ment of the DM induced gamma-ray signal over the one 
predicted based on the J factors calculated for the NFW 
profile and assuming the value of the velocity averaged 
cross section of (crv}th — 3 x 10~ 26 cm 3 s _1 . This value 
of the cross section leads to the correct relic abundance 
of DM in the simplest thermal decoupling scenarios and 
depends only logarithmically on the DM mass. We note 
here that large boost factors could originate both from 
an enhancement in the cross section or from an enhanced 
DM density. Non-perturbative electro-weak resonant ef- 
fects could naturally become relevant for heavy DM can- 
didates and can enhance annihilation cross sections, how- 
ever only for model-dependent narrow ranges of the DM 
mass, [76j . We will mention some other possibilities in 
the following section. 

The contraction of the dark matter density profile 
around the black hole residing in the center of the Galaxy 
could result in a so-called dark matter 'spike', a region 
of enhanced density. It could form due to the adiabatic 
growth of the black hole, due to scattering of dark mat- 
ter with the dense stellar environment around the black 
hole or due to baryonic infall [77M79| (however events 
such as binary black hole mergers or off-center super- 
massive black hole formation could results in the oppo- 
site effect.) In the simple adiabatic contraction model the 
slope of the spike is derived in terms of the inner slope 
7 as 7 S p = 2 + 1/(4 — 7) an( i nas an effect at distances 
below few parsecs from the black hole. Any of the above- 
mentioned mechanisms could provide an enhancement of 
the dark matter annihilation signal by a factor of 100 to 
1,000 in the innermost region of the Galaxy (for compar- 
ison of J-factors in such scenarios we refer the reader to 



B. Spectrum 

The spectrum of multi-TeV DM particles has to be 
treated with care as in collisions with center-of mass en- 
ergies significantly higher than the weak scale the radia- 
tion of electrowcak gauge bosons from the products of the 
DM annihilation/decay becomes possible. The so-called 
electroweak radiative corrections, neglected until recently 
[53|-[58|, are expected to have a sizable impact on the 
energy spectra and branching ratios of SM particles orig- 
inated from the annihilation/decay of DM particles with 
mass larger than the electroweak scale. 

Under the assumption that tree-level cross sections for 
processes leading to two final states are dominant over 
those corresponding to processes with three final states 
including a gauge boson it is then generically possible to 
compute such corrections in a model independent way 
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[54l| . These calculations have been implemented in the 
PPPC4DMID code [8l| 5 which we use to to derive the spec- 
tra of gammas and electrons from DM annihilations. 

A specific consequence of radiative corrections, rele- 
vant to our case, is that the total energy is distributed 
over a large number and variety of low energy particles 
amplifying the spectrum at lower energies (E Tn x ). 
We therefore do not expect them to impact our fits of 
gamma rays and electrons which are measured at high 
energies, close to the best-fit particle mass. However, 
positrons are measured at lower energies, where EW cor- 
rections affect the spectrum by up to an order of magni- 
tude dnn. 

Nevertheless, it should be mentioned that in some spe- 
cific cases radiative corrections can lead to large, model- 
dependent effects on the spectrum at high energies, with 
hard contributions near the DM mass. This can hap- 
pen for the emission of photons radiated from virtual 
charged particles and arc particular important for Ma- 
jorana DM or 'wino-like' scenarios [82l48^ |. These ef- 
fects lead to the appearance of a sharp line-like signa- 
ture, which could significantly increase the sensitivity 
of the DM searches [H, |86j]. However, the inclusion of 
these channels does not significantly improve the qual- 
ity of our gamma-ray fits (as also found in Re- 
cent works [H, [H, generalized the three-body final 
states calculations to also include the emission of EW 
bosons. These however affect the gamma-ray spectrum 
in the lower x = E^/moM ^0.3 region and are typically 
rather flat in energy [56] in this range. As the exact spec- 
tra and amplitude of these effects arc model-dependent 
we do not consider them here. 

In [87| it was shown that for the 'wino-like' DM candi- 
dates which are expected to mostly produce WW states 
via inclusion of photon emissions from the initial states 
actually can make comparable annihilations to three- 
body final states including fermions ffV (depending on 
the choice of couplings and exchange particle mass) . This 
is interesting in the view of our analysis here as, as we will 
see, the gamma ray spectrum is best fit with a combina- 
tion of WW and harder tt channels. This issue merits a 
dedicated investigation which is beyond the scope of this 
work. 

In our analysis, the sum runs over , t + t~ , bb and 

tt annihilation channels. As the DM annihilations to the 
bb channel result in a gamma-ray spectrum which is very 
similar to the spectrum of tt and W + W~ we will use it 
to mimic all non-leptonic annihilation channels. 



5 This calculation assumes that final states are massless, which is 
in general a good approximation given the high center-of-mass 
energies that we deal with here. 



IV. FITTING GAMMA RAYS FROM THE 
GALACTIC CENTER 

In this work we use the spectral analysis of Fermi LAT 
data from [39[ and HESS data from [22| ■ We study them 
in terms of an astrophysical component and a dark mat- 
ter particle annihilation spectrum, to a combination of 
representative channels (fi + fi~, t + t~, bb). Specifically, 
under the assumption that a power-law spectrum pro- 
vides a reasonable fit to point sources (see, e.g., the anal- 
ysis of the 1FGL catalog, |46|). our model for the joint HE 
and VHE spectra of the central source (1FGL J1745.6- 
2900) consists of the combination of a power-law and two 
additional contributions. 

As a possible alternative, we consider as well a 'log- 
parabola' spectrum (shown to better reproduce the spec- 
trum of AGNs in the 2FGL) which adds two parameters 
while decreasing smoothly at high energy. This will be 
given by 

N/dE = K(E/E o )- a - 0l ° 9(E/Eo) , (4) 

with parameters K, a (spectral slope at Eq) and the 
curvature (3, while Eq is an arbitrary reference energy. 

In analyzing the data, we consider three different sets 
of data. The first (A) consists of the Fermi LAT data 
excluding the last data point, assumed to be an outlier, 
and the HESS data, excluding the last four data points. 
These omissions are justified by their large error-bars and 
the potential difficulty in distinguishing the signal from 
the diffuse emission from the larger ROI in the more un- 
certain high energy region. In the second data set (B), 
we consider just the HESS data, still excluding the last 
four data points. The third case (C) is given instead by 
the full HESS and Fermi LAT data sets. 

The choice of these different combinations is motivated 
by our specific interest in multi-TeV cutoff apparent in 
the HESS data. Since the best-fit model would, to a 
large extent be determined by the Fermi LAT data which 
present much smaller errors, a comparison of the results 
from sets A and B is particularly relevant. In what fol- 
lows we will be mostly considering set A, unless otherwise 
specified. 

As mentioned above, we fit different functional forms 
to mimic possible astrophysical contributions in the GeV 
energy range and we add a DM spectrum to model the 
VHE range. In each case, we minimize the \ 2 over the 
mass of the dark matter candidate as well as the normal- 
ization of the DM flux, and the parameters of the addi- 
tional function(s). We also marginalize over the branch- 
ing ratios to , t + t~ and bb (with this last spectrum 
mimicking all non-leptonic annihilation channels). 

We start by fitting a power-law of the form Ae~ T 
alone (case G.l) and a combination of a power law 
with the spectrum produced by annihilating dark matter 
Ae~ T + <p x (E) (case G.2) to the data set A, see Fig. [TJ 
top left panel. These two cases clearly cannot accurately 
reproduce the data sets, and we quote values of x 2 m 
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case 


datasets 


r 


A xl(T 8 


GeV range fit 


DM mass [TeV] 


(Br T+T -; Br b - b ) 


d( 2 ) 
^NFW 


X 2 /d.o.f. 


G.l 


Fernii+HESS 


2.48 


6.45 


- 


- 


- 


- 


5.48 


G.2 


Fermi+RESS 


2.51 


6.63 


- 


12.7 


(1.0;0.0) 


390 


3.34 


G.3 


Fermi+RESS 


- 


0.25 


( E /4.2)- 2 ' 54 -°- 071 °s( B / 4 - 2 ) 


15.5 


(0.45;0.55) 


1645 


1.4 


G.4 


Fermi+RESS 


2.48 


4.8 


E -1.5 e -B/3.4 


12.2 


(0.8;0.2) 


507 


1.0 


G.5 


Fermi+HESS 


2.55 


4.4 


m x = 28, B F = 62 


18.2 


(0.48;0.52) 


1400 


0.9 


G.6 


HESS 


2.51 


3.21 




18.2 


(0.53;0.47) 


1075 


0.81 


G.7 


Ferrm+HESS(mll) 


2.59 


3.54 


£-1.44 e -E/3.8 


14.9 


(0.47;0.53) 


1322 


1.42 


G.8 


Fermi+HESS(full) 


2.51 


3.98 


28 


18.2 


(0.42;0.58) 


1685 


1.09 



TABLE I: Best-fit parameters and \ 2 f° r different models of the galactic center gamma-ray spectra. 



Table U just for reference. In the following we try to im- 
prove the low energy fit in order to gauge the impact on 
the DM fits in the high energy range. 

In the first place, we model the emission in the GeV 
region with an addition of a power law with an expo- 
nential cut-off to mimic the potential contribution of an 
unresolved pulsar population in this region. This will be 
given by 



E 



puis 



E 



_ 7 -E/E a 



(5) 



In this case we find that the best-fit parameters are 
E C ut = 3.8 GeV, 7 = 1.6, intriguingly similar to the mean 
values found in the Fermi LAT pulsar catalog 44|. We 
do not attempt to make strong statements about the ori- 
gin of this feature, as it might be challenging to localize a 
significant pulsar population in a tiny region of about 30 
pc 6 (pulsar kicks of 100 km/s allow them to move 10 pc 
over 0.1 Myr time, assuming they are not gravitationally 
bound), but we note that that such a spectrum provides 
a good fit in this case. An HST-based measurement per- 
formed in [HI suggests a high star formation rate in this 
region rich with young (few Myrs) stellar clusters which 
would naturally produce a numerous pulsar population. 
In addition older clusters could spiral into the central 
parsec region drawn by dynamical friction (89j (see also 
|45| for an interpretation of the extended source Sgr A* by 
millisecond pulsars). We denote the model that includes 
this exponential cutoff to the power-law and the annihi- 
lation spectrum of heavy dark matter fit to data set A 
as case G.4 in Table U (see Fig. [TJ bottom left panel). 
In addition, and following the Fermi LAT catalog 2FGL, 
we model the GeV emission by the log-parabola fit, case 
G.3. Wc find the best-fit values to be a = 2.5, /3 = 0.074 
and E = 4.2 GeV. This case is shown in the top right 
panel of Fig. [1] . 

Lastly, we model the low energy spectrum with a power 
law and an additional DM spectrum, which we again 
marginalize over all the DM masses and channels, case 
G.5 (Fig. [I] bottom left panel). In this way, we allow the 



most freedom to the fit, and this check is meant to test 
whether there is a significant impact on the heavy DM 
derived regions, by allowing for a large freedom in the fits 
to the low energy data. The \ 2 values for these choices 
are comparable as well as are the inferred properties of 
the heavy DM. 

In case G.8 wc consider the full Fermi and HESS 
datasets C (Fig. [TJ bottom right panel) and check that 
the fit parameters are not too different from the ones 
derived in the analysis of dataset A. 

As an additional check we consider as well an analysis 
limited to the HESS data (data set B), defined by a power 
law and a heavy dark matter annihilation spectrum (case 
G.6 in Table IJ). As mentioned before, a power law with 
an exponential cutoff instead of the superposition of a 
power law and a heavy dark matter contribution provides 
a comparably good fit to the HESS data with x 2 /d.o.f. = 
0.78. 

In the top panels of Fig. [5]wc plot the regions in the 
parameter space m x -Bp, as defined in Section [ill Al for 
the case G.4, marginalized over all DM annihilation chan- 
nels. These results are very similar to the results of case 
G.5. If we fit only the HESS data, the region gets some- 
what enlarged towards smaller boost factor values. In the 
lower panels we plot instead the allowed regions of param- 
eter space considering the DM mass versus the branching 



ratio to t + t , that is m x -Br T 



We assume the re- 



corresponding to the PSF of Fermi LAT at 3 GeV which is > 0.2 



maining annihilation products to proceed through softer 
bb-like channels. In both figures, light-blue and blue re- 
gions correspond respectively to 95% and 99.999% C.L. 
regions. 



V. LOCAL POPULATION OF ELECTRONS 
AND POSITRONS AND DM SPECTRAL FITS 

The standard established mechanism for producing 
cosmic ray (CR) positrons is by means of secondary pro- 
duction in the interactions of CR nuclei with interstellar 
gas. Their measurement was therefore recognized as a 
way to test models of the production and propagation 
of cosmic rays in the Galaxy, as well as models of new 
physics, in particular the annihilations or decays of dark 
matter 1901 . 
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FIG. 1: Fits to the gamma-ray spectrum of the source 1FGL J1745. 6-2900 using Fermi LAT and HESS data. Top, left panel: 
single power law, case G.l {dashed, blue curve) and a power law in combination with a heavy dark matter annihilation spectrum, 
case G.2 {solid, black curve). Top, right panel: log-parabola spectrum, modeling an average pulsar, and a heavy dark matter 
annihilation spectrum, case G.3 {solid curve). Bottom, left panel: a power- law and two dark annihilation spectra, case G.5 
{solid, black curve) and a power law with an exponential cutoff, modeling the pulsar contribution, case G.4 {dashed, blue curve). 
Bottom, right panel: same the bottom left panel but for data set C, cases G.7 {dashed, blue curve) and G.8 {solid, black curve). 
In all panels , the dark matter annihilation spectra, the power laws and the exponential cutoff function are plotted separately 
as dotted lines. See Table U for details. 
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FIG. 2: Top panels: allowed regions of parameter space m x -BF (mass of dark matter particle - boost factor) that fit Fermi 
and HESS galactic center gamma-ray data, case G.4 {top left panel) and that fit HESS data alone, case G.6 (top right panel) 
(the values of \ 2 f° r given values of m x and BF are found by minimizing over a combination of t + t~ and bb dark annihilation 
channels). Bottom panels: Allowed regions of parameter space m x -Br T + T - (mass of dark matter particle - branching ratio to 
t + t~) that fit Fermi galactic center gamma-ray data and HESS galactic center gamma- ray data, case G.4 (bottom left panel) 
and the HESS data alone, case G.6 (bottom right panel). In all cases, light-blue and blue correspond respectively to 95% and 
99.999% C.L. regions. 



This approach proved to be fruitful as the PAMELA 
satellite reported a surprising rise of the positron fraction 
e + /(e + + e - ) above the expected declining astrophysi- 
cal background from ~ 10 GeV up to at least 100 GeV 
[9l| (these data confirmed the results from the earlier 



HEAT balloon experiment [92[ and AMS test- flight [§3j]). 
Secondary production mechanism results in a positron 
fraction that decreases with energy, once diffusion effects 
are taken into account [94| . The rising positron fraction, 
also recently confirmed by Fermi LAT measurement, [95| . 
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therefore points to an existence of a source of primary 
positrons. 

At the same time, PAMELA measured anti-proton 
fraction, p/(p + p), consistent with expectations based 
on astrophysical secondaries, up to the maximal probed 
energy of about 200 GeV [9g|, in agreement with earlier 
data. 

Consistently with the positron rise, the measurement 
of the total electron (e + + e~) spectrum by Fermi LAT 
[97| showed a hard E~ 2 power law, opening the possi- 
bility that an additional hard spectral component con- 
tributes on top of a softer astrophysical spectrum. This 
measurement was later confirmed by PAMELA [98[ and 
updated by Fermi LAT [95|. The HESS telescope also 
reports the measurement of the e + + e~ energy spectrum 
above energies of 600 GeV, showing a power law spec- 
trum in agreement with the one from Fermi LAT and 
measuring a cut-off at around 3 TeV, [§9[. This spec- 
trum was subsequently c onfir med by the measurement 
by the MAGIC telescope [lOOj . 

There are several candidates for astrophysical sources 
which may supply this extra-yield of cosmic positrons: 
i) pulsars could provide sizable contributions to the 
positron flux from pair conversio ns of 7 ra ys in their mag- 
netic fields (for recent work see (lQlUlO^ h ii) spallation 
processes of cosmic rays i nside SNR s with the gas, dur- 
ing the acceleration stage [l05Hl07| or iii) a more refined 
treat ment of a spatial distribution of sources and gas 
jl08| . In analogy, dark matter annihilation (or decay) can 
also provide a needed source of positrons, provided that 
its annihilation cross section is enhanced with respect 
to the value expected for WIMPs and that dark mat- 
ter is either heavy, or has suppressed decay to hadronic 
states, ot herw i se th e antiproton measurements would be 
exceeded [l09l fTToj. 

In the comprehensive work of the authors mod- 

eled electron/positron contribution of pulsars and local 
supcrnovac, and showed that, while one cannot predict 
the local electron and positron fluxes (as the overall spec- 
trum is mostly set by a hierarchy in the local sources), 
it is still possible to explain current measurements by 
means of reasonable parametcrizations of the source and 
propagation, without requiring exotic contributions. 

In early calculations of the dark matter elec- 
tron/positron contribution to the locally measured spec- 
tra, energy losses were ca lculated assuming the Thomp- 
son regime [l IOL ll 12Ml 19| . However, it has been noted 
that the Thomson approximation is no longer valid for 
energies at Earth above a few tens of GeV, for which a ful l 
relativistic description is consequently necessary, |l20j [. 
Here we employ a full relativistic Klein-Nishina calcula- 
tion as implemented in the GALPROP code (see also 1 1 2 ll ] ) . 
We also use an injection spectrum of dark matter anni- 
hilation products which takes into account electroweak 
bremstrahlung 8l|, important for dark matter and spe- 
cially for positron fraction at lower energies. 

In what follows, we set aside the question of the nor- 
malization of DM signals in the GC and local regions 



and ask the question whether there is a spectral overlap 
between the gamma-ray feature and the feature in the 
electron and positron data at TeV energies. 

The data on the total flux of electrons and positrons 
are available in the range 1 GeV through 4.5 TeV, with 
PAMELA covering the lower part from 1 through 625 
GeV [13, Fermi LAT in the range 20 GeV-1 TeV H3 
and HESS [H in the range 370 GeV - 4.5 TeV. To ac- 
count for systematic uncertainties to the statistical errors 
provided in the HESS dataset, we have added a system- 
atic uncertainty from the HESS band bin by bin. 

We will use the term 'electrons' to imply the total elec- 
tron plus positron flux which we describe by the sum of 
a power-law astrophysical component (modulated by the 
solar potential when required) plus a dark matter contri- 
bution. We note that a simple power-law is not expected 
to necessarily be the best fit for the astrophysical electron 
spectrum in the high energy range, where the stochastic 
nature of nearby sources needs to be taken into account 
but we resort to this parametrization for simplicity. In 
addition, to circumvent this problem while exploring the 
position of the energy (or DM mass) cut-off, we focus 
only on the HESS data. 

We have found the best fits for the combination of 
the PAMELA, Fermi LAT and complete HESS datasets 
modeling the attenuation of the flux at energies below 10 
GeV due to solar modulation by rescaling the spectrum 
[l22l Il23| . A first approach consists in keeping the solar 
modulation potential as a free parameter (case E.l), 
while a second approach assumes a fixed value <p = 600 
MV (case E.2). To minimize the effect of solar modula- 
tion, we have also considered a power law and dark mat- 
ter fit for the Fermi LAT and the complete HESS datasets 
alone, thereby restricting ourselves to the range above 20 
GeV (case E.3) but including a fit to the positrons data 
from Fermi LAT. Finally, case E.4 corresponds to fitting 
just the HESS datasets with dark matter. In Table |H] 
we list the parameters values corresponding to the best 
fits for all mentioned cases. The results for the cases 
mentioned above are shown in Fig. [3] 



Fig. 0] shows the regions of parameters space m x - 
Br T + T - (dark matter mass - branching ratio to t + t~) 
that fit PAMELA, Fermi LAT and HESS electrons and 
positrons data with varying solar potential <fi, case E.l 
(left panel) and e + + e~ data with, in addition, the Fermi 
LAT e + data, case E.3 (right panel). Light-orange and 
red correspond respectively to the 95% and 99.999% C.L. 
regions. 

In Fig. we plot the allowed regions in the same 
parameter space m x -Br T + T - from the HESS GC gamma- 
rays analysis, case G.6, and HESS e + + e~ data, case E.4. 

Finally, Fig. [5] shows the regions in the parameters 
space m x -boost factor for the e + + e~ data corresponding 
to case E.3. The boost factors required by the DM fits 
in this case are an order of magnitude higher than the 
ones implied by the GC gamma-ray data. This poses 
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case 


datasets 


mass[TeV] 


(Br T+T -; Br b - b ; Sr H+H -) 


10~ 3 x Bnfw 


r 


A 


4> [MV] 


X 2 /d.o.f. 


E.l 


PAMELA+Fermi+HESS 


3 


(0.9997;0;0) 


12.8 


3.72 


1943 


986 


0.3 


E.2 


PAMELA+Ferrm+HESS 


3 


(0;0;0.9996) 


1.38 


3.33 


595 


600 


1.3 


E.3 


Fermi +HESS (+ Fermi e+) 


3 


(0.9997;0;0) 


8.7 


3.32 


400 




0.05 


E.4 
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94.9 








0.4 



TABLE II: Fit to the electron and positron spectra 
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FIG. 3: Left panel: best fits for case E.l (solid curve) and E.2 (dashed curve). Right panel: best fits for case E.3 (solid curve) 
and E.4 (dashed curve). In both panels, individual components of the fits (power law and the dark matter spectrum) are shown 
by the dotted lines. The best-fit values for the parameters are given in Table ITT1 



an additional difficulty in relating these two phenomena 
and at the same time points to the fact that a potential 
signal in the GC regions would not be in a conflict with 
the local measurement of electrons/positrons. 

In order to alleviate the problem of high boost factors, 
an explanation of the electron flux in terms of dark mat- 
ter annih ilating in a nearby clump has been proposed in 
jl24ll25j . The normalization of the flux then depends on 
the size and the density profile of the clump as well as on 
the distance to the clump. The farther away we place the 
clump, the bigger and/or denser we have to make it to 
achieve the same required boost factor. The effect on the 
spectral shape, however, is not degenerate. In fact, the 
high-energy cutoff in the electron spectrum is inversely 
pro porti onal to the distance to the source of electrons 
[fjfjl Il26l | . In such a way a spectrum of electrons which 
initially extended to higher energies would appear to us 
as being cut off to lower energies and it is conceivable 
that an explanation of electron flux in terms dark matter 
annihilating in a nearby clump could be more spectrally 
consistent with the explanation of gamma-ray data from 
the galactic center. 

However, large clumps needed are very unlikely at the 



required dista nce t o the Solar neighborhood |l27l | and it 
was shown in [114| that the average contribution of the 
clumps is expected to be negligible compared to that of 
the smooth dark matter component. Despite these con- 
siderable difficulties, a study of electrons due to clumpy 
DM in the light of models we consider here could be of 
potential interest. 



VI. CONSTRAINTS ON ANNIHILATION 
CROSS-SECTION FOR HEAVY DM 
CANDIDATES 

In this section, we briefly compare our results with cur- 
rent constraints, focusing on heavy ( > 1 TeV) DM candi- 
dates (for a recent review on indirects constraints on DM 
annihilation cross section see (fjfjj.) We should underline 
that the analysis in this work relies mainly on spectral 
features, whose normalization constrain the product of 
cross section and dark matter density. The significant 
uncertainty on the density in this particular region re- 
sults in turn in a large uncertainty on our predictions for 
the annihilation cross section ((<rv)) so that comparisons 
with the limits on (cry) listed below should therefore be 
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FIG. 4: Allowed regions of parameters space m x -Br T + T - (DM mass - branching ratio to t + t~) from the analysis of PAMELA, 
Fermi LAT and HESS e + + e~ data, case E.l (left panel) and Fermi LAT and HESS e + + e~ data with, in addition, the Fermi 
LAT e + data, case E.3 (right panel). Orange and red correspond respectively to the 95% and 99.999% C.L. regions. 
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FIG. 5: Allowed regions of parameter space m x -Br T + T - (DM 
mass - branching ratio to t + t~) from the fit of the HESS GC 
gamma-ray data, case G.6 and the HESS e + + e~ data, case 
E.4. Light-blue and blue correspond to the 95% and 99.999% 
C.L. regions from the gamma-ray data and orange and red 
correspond instead to the 95% and 99.999% C.L. regions from 
the electron data. 



taken only approximately. 

• HESS GChalo analysis [H| (see also PH). HESS 
observations of the Galactic Center present the 
strongest constraints on WIMP dark matter anni- 
hilation into Standard Model particles for > 1 TeV, 
given a dark matter Milky Way density described 
by a non-adiabatically-contr actc d NFW or Einasto 
profile. DM signal search in jl28l | is performed with 
a projected galactocentric distance of 45 pc - 150 
pc (corresponding to an angular distance of 0.3° - 
1.0°), excluding the Galactic plane. Note, however, 
that for an isothermal halo profile, the signal would 
be completely subtracted in this analysis (as back- 
grounds are modeled according to the signal in the 
sorrounding region) and therefore if the DM signal 
is flat at distances > 45 pc, this limit does not ap- 
ply. Limits are derived under the assumption that 
DM particles annihilate into quark- ant iquark pairs, 
and constraints DM annihilation cross sections of 
the order 2 x 10~ 24 cm 3 s -2 , or boost factors of 140 
in our notation (when we renormalize to the same 
local density of DM) at 10 TeV. 



GC ridge [1301 ]. A band of diffuse emission along 
the GC ridge, was used to de rive upper limits on a 
DM induced signal The limits are some- 

what worse than in the analysis mentioned above, 
and are of the order of 10 -23 cm 3 s~ 2 , or boost fac- 
tors of 700 in our notation at 10 TeV, for the an- 
nihilation channel to t + t~ . 
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PAMELA antiproton data [96| were used to set lim- 
its on DM annihilat ions to h adronic final states 
[lolllioj]. Recently [HI Eli considered the case 
in which the dark matter is constituted by Majo- 
rana particles and studied the annihilation process 
into two fermions and one gauge boson which, un- 
der some circumstances, can have a non-negligible 
or even a larger cross section than the lowest order 
annihilation process into two fermions as we dis- 
cussed in section IIIIBI They performed a general 
analysis of this scenario, calculating the annihila- 
tion cross section of the three body final state pro- 
cesses, XX ~ ► ffV when the dark matter particle is 
a SU(2)l singlet or doublet, where / is a lepton or 
a quark, and V is a photon, a weak gauge boson or 
a gluon. This framework is relevant in the context 
of this analysis as heavier DM is expected to pro- 
duce an additional gauge boson more readily.They 
set a limit on the value of the boost factor at ~ 10 3 
at 5 TeV, for the case of singlet DM coupled to 
quarks (limits in the case of a coupling to leptons 
are weaker) . 
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boost factor over NFW for r r and 66 channels vs m v , case E. 3 



Neutrino observ ations (I CE CUBE [136j and Su- 
perKamiokande jl33lll37| b Such observations have 
the advantage that higher energy neutrinos coming 
from heavier DM annihilations also have a higher 
cross section for detection and thus the loss in rate 
is partially compensated by the increase in sensitiv- 
ity. Neutrino constraints become therefore some- 
what competitive with 7-ray bounds at large DM 
masses, reaching cross sections of ~ 10~ 22 cm 3 s -2 
at 10 TeV. 

Multi-wavelength analysis. A systematic study of 
the multi-wavelength signal induced by WIMP an- 
nihilations at the GC was performed i n IHTI (se e 
also jl38| for X-ray constraints, and [l39l . fl40j ] . 
for synchrotron constraints). The authors shown 
that the luminosity of the WIMP source and ex- 
perimental sensitivity at different frequencies is at 
a comparable level, depending rather weakly on 
WIMP mass and annihilation channels. For 10 TeV 
DM mass and t + t~ annihilation channels the most 
stringent constraints come from the near-infrared 
VLT data in the Ks band (2.16 /im) based on 
the observation within the central 10-20 milliarcsec 
(<0.01 pc). For the assumed adiabatically con- 
tracted DM profile (enhanced by about a factor of 
100 with respect to the NFW profile in that region) 
limits are at the level of (av) < 10~ 26 cm _3 s -1 , 
roughly consistent with the ones derived in this 
work. These limits however depend on the details 
of the DM density profile as well as on the energy 
loss and propagation parameters of electrons pro- 
duced in DM annihilations in the close vicinity of 
the black hole. 
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Unitarity arguments. In |141| it was shown that 
unitarity arguments imply interesting constraints 



FIG. 6: The regions of parameter space m x -BF (mass of dark 
matter particle - boost factor) that fit Fermi and HESS elec- 
tron data, case E.3 (the values of \ 2 f° r given values of m x 
and BF are found by minimizing over a combination of r + r _ 
and bb dark annihilation channels). Faded orange and red 
correspond to 95% and 99.999% C.L. regions. 



on strongly annihilating dark matter constrain- 
ing annihilation cross sections to be weaker than 
- 10~ 21 cm 3 s~ 2 at 10 TeV. 



VII. SUMMARY 

In this work, we fit the gamma ray spectrum of the 
central Galactic source over five decades in energy from 
the joint Fermi LAT and HESS data by making several 
choices to model the astrophysical emission and adding 
a heavy DM contribution. In particular we model the 
joint spectra with a combination of a power law, a spec- 
trum of unresolved pulsars and a DM contribution. Good 
fits are achieved with a power law with a slope index of 
~ 2.5 — 2.6, an exponential cut-off with index similar to 
those of the observed pulsar population and > 10 TeV 
DM decaying to a mixture of b b and harder t + t~ chan- 
nels. We show that, alternatively, the GC gamma-ray 
spectrum can be well fit with a log-parabola and heavy 
DM spectrum alone. The best fit in this case is real- 
ized by a somewhat heavier DM annihilating with higher 
branching ratios to channels resulting in softer f b b- like 
spectra. We also performed a fit to the HESS data alone, 
and conclude that the DM parameter space is less con- 
strained in this case, allowing for a wider mass range and 
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combinations of hard and soft channels. 

The best fit DM models generally require boost fac- 
tors of the order of 100-300 (for a DM mass of 10 TeV), 
when compared to the values expected based on the NFW 
profile and a thermal annihilation cross-section. These 
values are consistent or in a mild tension with the con- 
straints derived analyzing the HES S data in the region 
45-100 pc away from the GC [l28j |. However, the DM 
profiles in the vicinity of the super massive black hole 
could be slightly enhanced with respect to their shape 
at distances larger than 50 pc away, where GC halo con- 
straints were derived, thereby loosening such constraints. 

In addition, by leaving a large freedom in our param- 
eters, we explore whether there might be an overlap in 
the DM parameter space obtained when considering the 
Galactic center gamma ray and local electron/positron 
data. We studied the DM mass-branching ratio to r al- 



lowed parameter space, m x -Br T + T - , for the gamma ray 
and electrons and positrons data, based on spectral fea- 
tures and find that the 99.999% C.L. regions overlap only 
marginally. 
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